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Slow relaxation occurs in many physical and biological systems. 'Creep' is an example from 
everyday life: when stretching a rubber band, for example, the recovery to its equilibrium length 
is not, as one might think, exponential: the relaxation is slow, in many cases logarithmic, and can 
still be observed after many hours. The form of the relaxation also depends on the duration of 
the stretching, the 'waiting-time'. This ubiquitous phenomenon is called aging, and is abundant 
both in natural and technological applications. Here, we suggest a general mechanism for slow 
relaxations and aging, which predicts logarithmic relaxations, and a particular aging dependence on 
the waiting-time. We demonstrate the generality of the approach by comparing our predictions to 
experimental data on a diverse range of physical phenomena, from conductance in granular metals, 
to disordered insulators, and dirty semiconductors, to the low temperature dielectric properties of 
glasses. 

Keywords: Aging — Glasses — Slow relaxations 



Physicists often take for granted that systems relax 
exponentially. Indeed, when a capacitor discharges, it 
will discharge exponentially, with a rate independent of 
the time it has been charged for. However, the relax- 
ation of many systems in nature is far from exponential, 
as was noticed already in the 19th century by Weber 
In many cases, the relaxation is logarithmic: such 
relaxations have been experimentally observed in the de- 
cay of current in superconductors [^l, current relaxation 
in MOSFET devices mechanical relaxation of plant 
roots [H, volume relaxation of crumpling paper (Hj and 
frictional strength to name but a few. Fig. [T] shows 
experimental data for electron glasses and of crumpling 
a thin sheet, that are governed by extremely different 
physical processes, yet they display identical relaxation 
behavior, which is logarithmic over a strikingly broad 
time window. 

In these systems, in contrast to the capacitor example, 
the relaxation does depend on the time the system has 
been perturbed for ~ in the scientific jargon, this is re- 
ferred to as 'aging'. In fact, slow relaxations and aging 
are amongst the most distinct features of glasses, whose 
understanding presents an important problem in contem- 
porary condensed matter physics. Much experimental 
and theoretical attention has been devoted to aging in 
the past decades, in a variety of fields, such as spin-glass 
0] , colloids Q , vortices in superconductors Q and many 
others 

Here we study a generic model for aging, and discuss 
several mechanisms yielding a broad distribution of relax- 
ation rates. We demonstrate the generality of the model 
on four different experimental systems, measuring the 
dependence of the relaxation both on time t and on the 
'waiting-time' tw , during which an external perturbation 
has been applied. We show how the following form of 



relaxations transpires: 



Sit, t^) (x\og{l +t^/t) = 



log(i„,/t) for i <C 



t^/t 



for t > 



(1) 



where S is the physical observable. This means that the 
initial relaxation, at times short compared with t^, is 
logarithmic, while at long times compared to t^, it falls 
off as the inverse of the time - a power-law decay, much 
slower than exponential or stretched exponential decay 
p. 

Fig. [2]shows the excellent agreement between this pre- 
diction and experimental results for four different sys- 
tems, measuring various physical observables: conduc- 
tance relaxation in the electron glasses indium oxide and 
granular aluminum [l^ - l2]| . relaxation of the dielectric 



constant in the plastic mylar |22h24|| and conductance 
relaxation in room temperature porous silicon 25, [26| . 
The experiments also markedly differ in the involved 
timescales. Details of the experiments are given in Ta- 
ble 1. 

In the following, we describe the aging protocol used 
in the experiments, and introduce a model that predicts 
Eq. 11]). We explain how one can understand the slow 
relaxations in terms of an underlying distribution of re- 
laxation rates of a particular form: 



P(A) ~ 1/A, 



(2) 



which we shall show can emerge due to various, physi- 
cally distinct, mechanisms: thermal activation, quantum 
mechanical tunneling, or through a third mechanism re- 
lying on a multiplicative process. We then proceed to 
describe the connection to the ubiquitous 1// noise en- 
countered in many systems, as well as the possible role of 
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FIG. 1: Experimental results showing a logarithmic relaxation in the electron glass indium oxide, where conductance is 
measured, and in a system of crumpling mylar, where the height is measured, after a sudden change in the experimental 
conditions. As seen in the graph, the logarithmic change in the physical observable can be measured from times of order of 
seconds or less to several days [a, Similar logarithmic relaxations, observed over many decades in time, occur in numerous 
physical systems, ranging from currents in superconductors to frictional systems. Data courtesy of Z. Ovadyahu and S. Nagel. 



TABLE I: Details of experiments 

System Measured variable Units (s) Taken from 

Aluminum Conductance (a) 0.025a/a 1200960, 360 Grenet et al., Ref. [19] (2007). 

Mylar Dielectric constant (e) W'^Se/e 18000 Ludwig et al., Ref. [22] (2003). 

Indium Oxide Conductance (a) 0.015a /a 20, 180, 1620 Ovadyahu, Ref. [15] (2006). 

Porous Silicon Conductance (a) 0.02Sa/a 110 S. Borini, private comm. (2010). 
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the distribution described in Eq. (f5| in other intriguing 
phenomena such as Bcnford's law [27| . 



THE EXPERIMENTAL PROTOCOL 

The aging protocol is illustrated in Fig. [3l Its first 
step consists of letting the system attempt to equilibrate 
for a relatively long time (typically of the order of hours 
or days). Next, one perturbs the system, in a way which 
depends on the experimental system: for indium oxide 
and aluminum this is done by changing the voltage of a 
capacitivcly coupled gate, for the mylar sample it is done 
by putting the system in a perpendicular electric field, 
while for porous silicon a large bias voltage is applied. 
The perturbation is now maintained for a time . After 
it is switched off, the physical observable is continuously 
monitored, as it relaxes. The longer is, the slower the 
resulting relaxation. 



V, containing all the relevant information determining 
the physical observable. Perturbing the system weakly 
near its equilibrium, one can always linearize the equa- 
tions of motion, and obtain an equation = ASv, with 



Sv = v{t) — Vequiiibriumi and A a matrix, independent 
of time. The real part of the eigenvalues of A must be 
negative, in order for the equilibrium to be stable, and 
they have the physical meaning of relaxation rates: As 
is seen by solving the linear equation, each eigenmode 
relaxes exponentially to zero, independently of the other 
eigenmodes. In certain cases [ig, the form of A can be 
worked out explicitly. 

We shall now explain three different mechanisms which 
lead to an abundance of slowly relaxing modes of the 
system, described mathematically by Eq. Using our 
assumption that the modes contribute positively and uni- 
formly to the measured physical observable, the superpo- 
sition of these modes will yield (in a certain time window, 
the conditions of which we will discuss) the logarithmic 
relaxations described earlier. 



THE MODEL 

Having laid out a concrete experimental protocol we 
are now well positioned to describe a generic model which 
will yield the form of aging described by Eq. ([1]). The 
model will involve two ingredients: First, the understand- 
ing that a broad distribution of relaxation rates A occurs 
whose logarithm is approximately uniformly distributed 
over some broad range, as described by Eq. The 
second ingredient involves understanding what this relax- 
ation rate distribution implies on the relaxations, within 
the aging protocol. 

The model assumes that the measured physical observ- 
able (conductance, dielectric constant etc.) is affected by 
an ensemble of modes, which are independent and con- 
tribute to the observable in a definite way (for example, a 
relaxation of any of the modes will cause the conduction 
to decrease). For each particular system under study, 
understanding the microscopic source of these modes is 
a quite different sort of question one may ask, and is 
outside the scope discussed here |28|. In general, these 
assumptions would make sense for a physical observable 
which depends on the configuration of the whole system 
(e.g: conductance, volume), and not some local probe 
(e.g: current measured in an STM tip). Also, wc would 
always be assuming that a large number of these modes 
contribute, so that wc can take the continuous limit, and 
discuss probability distributions. Presumably, the model 
would fail for a sufficiently small sample (although, in 
the case of electron glasses the logarithmic relaxations 
were exp erimentally observed even for micron sized sam- 
ples [221). Each of the modes would relax exponentially 
to its equilibrium, but with a different relaxation rate. 
To support these assumptions, one may think of the sys- 
tem as being formally characterized by a state vector 



Thermal activation 

A diversity of physical processes are governed by ther- 
mal activation, which is perhaps the simplest physical 
mechanism which can give rise to Eq. (I2|), as has been 
known for long [3d - l32| . We should have in mind a rugged 
energy landscape characterizing the system, with many 
local minima. At a given time, we can denote by p the 
vector of probabilities for the system to reside in each of 
the minima. Clearly, for a system governed by stochastic 
dynamics, the probability vector would obey the same 
linear equation mentioned earlier, namely, ^ = Ap [i.e., 
wc have defined a Markov process). The relaxation 
modes in this case approximately correspond to cross- 
ing one of the energetic barriers connecting two of the 
local minima. Here, the rate A of a given process is given 
by the Arrhenius formula, namely, it is exponential in 
the energetic barrier U, namely A oc e^^^'^'^ , which the 
system has to cross in order to reduced its energy. We 
will associate each mode with one such transition (across 
an energetic barrier U) , and assume that the size of these 
barriers is distributed smoothly over a certain range of 
energies. We can now readily calculate the resulting re- 
laxation rate distribution: 



P(A) = 



\dX/dU\ 



T/A, 



(3) 



where we have taken P(f) as approximately constant. 
For isothermal processes, the temperature dependence 
entering the proportionality constant does not play a role 
in the aging behavior. From this formula we can also de- 
duce the smallest and largest rates the system supports 
(corresponding to the fastest and slowest times): these 
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FIG. 2: (a) Results of aging experiments for four different systems measuring different physical observables. Experimental 
parameters can be found in Table HI The x axis denotes the time (on logarithmic scale, spanning five decades in time, from 
seconds to days), and the y axis denotes the signal (with different units for each data set), (b) The x coordinate of each data 
set is scaled according to the known waiting-time tw, and the y coordinate is scaled such that the signal at t = tw is log(2) (for 
convenience). The data collapses onto a single curve, which is compared to the theoretical prediction of Eq. ([l]). The inset 
shows the same data has indeed linear dependence when the x axis is defined according to Eq. ([l|. 



5 



extremal barrier heights. Note that due to the expo- 
nential dependancc on U, even a small range of energy 
barriers can result in a broad range of relaxation rates. 
The above mechanism is essentially the same leading to 
1/f noise 33|, |3J], and is reminiscent of Bouchaud's trap 
model [32 1 . 

It should be emphasized that we have assumed here 
that the energy barriers vary sufficiently slowly, and 
therefore their variation within the energy interval 
[Umim Umax] Can be neglected. This simple picture turns 
out to be extremely successful when applied to recent 
experiments on porous silicon performed around room 
temperature. At such a high temperature, thermal ef- 
fects are expected to be dominant over quantum effects, 
and indeed the maximal timescale Amax was found experi- 
mentally to be very sensitive to temperature, as expected 
from the above formula (26j . 

However, it is experimentally found that in various sys - 
tems the relaxations are insensitive to temperature [l9[, 
which necessitates a different mechanism. Quantum Tun- 
neling (QM) is the second mechanism which yields Eq. 
@. 



Quantum mechanical tunneling 

Let us keep in mind the picture described earlier for 
the rugged energy landscape, but now assume that we 
are at low enough temperatures such that thermal ac- 
tivation across the barriers is prohibited. The system 
will be able to quantum mechanically tunnel through 
the barriers, paying a penalty which is typically expo- 
nentially suppressed with the distance. Thus, for this 
process as well, the rate A depends exponentially on a 
smoothly distributed variable, which in this instance is 
the distance: A ~ e~^^^^, with ^ the localization length 
of the wavefunctions, and r the s pat ial distance between 
the two points. In a recent work |35| . the distribution of 
relaxation rates was calculated for this case, taking into 
account the correlations that exist (the distances are not 
independent in this case). It was found that one still 
obtains the 1/A distribution, albeit with small but inter- 
esting logarithmic corrections. Here, the role played by 
the temperature T in Eq. ([3]) is played by the localization 
length ^. Related considerations for a varying height of 
the barrier through which the tunneling occurs are given 
inRef. 

In both cases, of thermal activation and of quantum 
mechanical tunneling, we would like the barrier distribu- 
tion to be broad in energy or real space, namely, it should 
be large compared to the energy kT or the localization 
length ^, respectively, in order to achieve a broad range 
of relaxation rates. 

So far we discussed two different natural ways which 



cesses is the key ingredient in obtaining Eq. ([5]). Both 
mechanisms, however, are inadequate to describe the log- 
arithmic relaxation in crumpling paper, for example [5[. 
We shall now present another mechanism, which does 
not rely on a variable being exponential, but rather, on 
the central-limit-theorem. This suggests the mechanism 
should be broadly applicable. We will show how the in- 
terplay of many random processes can under general con- 
ditions lead to a log-normal distribution, which is well 
approximated over a broad range by Eq. ([2]). 

Multiplicative processes 

In many physical examples, an observable depends on 
the product of many approximately independent vari- 
ables, which is referred to as a multiplicative process. 
Understanding the importance of su ch p rocesses in na- 
ture dates back (at least) to Shockley [STI, who discussed 
the connection of multiplicative processes to log-normal 
distribution which we shall also utilize here. See [ssf for 
a strongly related discussion in the context of 1/f noise. 

An example of a multiplicative process is the transmis- 
sion of a particle through a one-dimensional disordered 
wire: if we divide the wire into a large number of seg- 
ments, it is known that the average transmission is the 
product of the individual transmissions through each seg- 
ment [3^. Fig. 4 demonstrates pictorially another such 
example, namely, how for electron glasses relaxation can 
occur via the simultaneous tunneling of various electrons, 
which is also approximately a multiplicative process. If 
we assume that the relaxation rate A is a product of many 
independent variables Xi , we can readily calculate the re- 
sulting distribution of relaxation rates. Since A = Yii ^ij 
we have, using the central limit theorem: 



P (log(A)) ^ e 



-[log(A)-,.]/A^ 



(4) 



By changing back to the variable A, we find that it follows 
a log-normal distribution: 



P(A) 



,-[log{Ato)/A] = 



/A, 



(5) 



lead to it, namely, thermal activation |30l - l32l | and quan- 



tum tunneling 35[ . The exponential nature of these pro- 



where is a constant with the dimensions of time. 
Far from the tails of the distribution, namely, when 
|log(Ato)| <C A the distribution reduces to Eq. Re- 
markably, in the crumpling paper example, the distribu- 
tion of the lengths of the segments was measured directly, 
and shown to follow a log-normal distribution [i^l- The 
compatibility of log-normal distribution and logarithmic 
relaxations fits well with the theoretical framework we 
suggest. 

DERIVATION OF THE AGING FORMULA 

We shall now discuss the implications of this distribu- 
tion on aging experiments, showing it leads to aging of a 
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FIG. 3: Schematic description of the different stages of the aging protocoL At time t = — t^, a perturbation is apphed to the 
system, which is turned off at time t = 0. We will be interested in the form of the relaxation of the physical observable in the 
last stage. 



(a) (b) 




FIG. 4: Pictorial demonstration of different physical mechanisms leading to a broad distribution of relaxation rates, (a) The 
energy landscape of various complex systems, including glasses, contains many minima. The energetic barriers separating them 
are smoothly distributed over a certain range. In order for the system to relax its energy, it must cross these barriers by thermal 
activation or by quantum mechanical tunneling: both are exponential in the barrier, which lead to a rate distribution described 
by -P(A) ~ 1/A, as we discuss in detail, (b) Many particle transitions in an electron glass are an example of a multiplicative 
process: In many electronic configurations, moving any single electron in the system to one of the vacant sites will result in 
higher energy, and therefore these processes will not occur at low enough temperatures. However, changing the position of a 
larger number of electrons can result in a lower energy. The rates of this process can be approximately written as a product of 
the rates of the single particle processes involved, leading to a rate distribution described by a log-normal distribution. 
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particular form, called 'full' aging, in a particular limit. 
This was done in the context of relaxation in electron 
glasses in [2^. The derivation we will present, however, 
does not rely on any peculiar properties of this specific 
system, and as such will be broadly applicable for all 
physical systems which follow the 1/A distribution. 

Let us assume that the system supports i = 1 . . . N 
1 relaxation modes, with corresponding relaxation rates 
Xi , each of which contributes an amount Xq to the physi- 
cal observable {e.g.: conductance or dielectric constant). 
Before going to the more involved aging experiment, let 
us consider the case where we excite all of these modes 
by some uniform amplitude. The physical observable X 
measured a time t after the perturbation, would read 
X{t) = XoJ2iS~^^* : which in the continuous limit goes 
to: 



X{t) = Xo 



dXP{X)e 



(6) 



where P(A) is the distribution of relaxation rates. In the 
case where P(A) ~ 1/A, introduced earlier, with Amin and 
Amax the lower and upper cutoffs, we obtain the difference 
of two exponential integral function IS 41 j : 



X{t) — Xo[i?l(Aniini) — £'l(AmaxO]- 



(7) 



For the case where the involved times are much smaller 
than the reciprocal lower cutoff Amin, and much larger 
than the reciprocal upper cutoff, the equation reduces to 
a simpler form: 



X{t) = Xo[-7£; -log(Ami„i)]- 



(8) 



Thus, we expect a logarithmic relaxation, which is indeed 
experimentally observed in a large variety of systems, 
as discussed earlier. Fig. 1 shows such a logarithmic 
relaxation, measured in an indium oxide sample. 

Going on to the aging protocol, we shall assume that 
the perturbation is small enough such that the rates 
of the relaxing modes arc indifferent to it. Nonethe- 
less, upon the application of the perturbation, the fixed 
point to which the system attempts to relax to (which 
does not have to be the true equilibrium, but can be a 
metastable state), is different when the perturbation is 
applied. Therefore during the second stage of the exper- 
iment (see Fig. 4), the system relaxes towards the new 
metastable state, which means that the relaxation modes 
are excited with respect to the original metastable state: 
relaxation to the new metastable state implies excitation 
with respect to the old one. The closer we got to the new 
minimum, the further we are from the initial one. 

Let us illustrate this for the example of a single relax- 
ation mode: in this case the relaxation is exponential. 



and therefore at the moment when the perturbation is 
switched off the distance from the new metastable state 
is proportional to e"^*™. At this moment the distance 
from the original metastable state is 1 — e~'^*™ : indeed, 
for — Q nothing happens, while for t^^, — > 00 the largest 
possible excitation occurs. A time t later, the ampli- 
tude of the relaxation mode, which decays exponentially, 
would therefore be (1 — e~'*'*™)e~'^*. Generalizing this for 
the case of many relaxation modes, we find that a time t 
after the perturbation has been switched off, the physical 
observable is given by: 



X{t) = 



dAP(A)(l 



-At 



(9) 



This can be written, as before, in terms of exponen- 
tial integral function. Assuming that we are in the 
regime of intermediate asymptotics, where the experi- 
mental timescales are much larger than the reciprocal 
upper cutoff and much smaller than the reciprocal lower 
cutoff, we obtain the difference of two logarithms: 



X{t)/Xo = log[Ami„(i + t^)] - log[Ami„i], 



(10) 



leading to Eq. dT]). 

It should be noted that the regime where the experi- 
mental time is comparable to l/Amin can also be reached, 
and it was shown that Eq. Q accounts of the aging be- 
havior in porous silicon also in the case where significant 
deviations from the full aging regime were observed [2^ . 
In other words, the model predicts full aging only in the 
asymptotic regime, and can also account for the devia- 
tions from full aging. Related system dependent cutoffs 
were also discussed in the context of spin-glass [42 1. 



CONNECTION TO 1// NOISE 

The broad underlying distribution of relaxation rates 
which played a crucial role in determining the slow relax- 
ations, is also central for understanding low- frequency 
noise in many systems. A variety of physical, biologi- 
cal and financial models show a universal form of low- 
frequency noise 3^, 34, 3^, 4^, with a power-spectrum 



scaling as 1 //. This ubiquitous form of noise is deeply re- 
lated to the logarithmic relaxation which we study here, 
the underlying principle bein g a roughly uniform distri- 



bution of effective barriers. In [4J|, a relation between the 
two physical phenomena was made, based on a theory de- 
vised by Onsager nearly a century ago: the connection is 
made through the Onsager's regression principle, stating 
that the relaxation of a system close to its equilibrium is 
related to the spectrum of the fluctuations of the system 
around the equilibrium [45|. Each mode with a rate A, 
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generates a Lorentzian noise spectrum Hill 



A 



/2 + A2 ■ 



(11) 



Summing over many modes using the distribution of Eq. 
(21); yields 1// noise. Linking these two seemingly un- 
related universal behaviors, logarithmic relaxations and 
1// noise, seems to us both conceptually appealing and of 
practical importance. For one, it means that the different 
physical mechanisms wc suggested to yield the broad dis- 
tribution of relaxation rates, would also imply 1// noise. 



OTHER APPEARANCES OF THE 1/A 
DISTRIBUTION IN NATURE 



In fact, the P{X) ~ 1/A distribution which plays a 
pivotal role in determining the aging, logarithmic relax- 
ations and 1/f noise, appears also in completely different 
contexts. The Gutenberg-Richter law [461 , for example, 
states that the distribution of the magnitude of earth- 
quake is a power-law. The exponent is experimentally 
found, in various methods of analysis, to be close to one 
[4?! . Another striking example lies in the so-called 'first- 
digit problem'. 

Towards the end of the 19th century, the astronomer 
Simon Newcomb noticed, while looking at his logarithm 
table, that numbers starting with 1 are looked up far 
more often than higher digits 



48|. Half a century 



later, the physicist Frank Bcnford rediscovered the phe- 
nomenon, and asked himself the following question: what 
is the distribution of the leading digits of numbers en- 



countered in a certain scenario? j27l |. Remarkably, in 
hugely differing data sets such as those found in tax re- 
turns, tables of physical constants, birth rates and many 
others, the relative occupance of each digit follows a uni- 
versal (and nonuniform) distribution: P{d) = log(l + 
l/d), where d is the digit. This is known as Benfords 
law, and is robust enough to be used to detect frauds in 
tax returns It can be explained in a simple way, if 
we assume that the distribution of x, the variable mea- 
sured, follows approximately P{x) ^ 1/x, over a large 
window. Clearly, the probability that the first digit is 1 
is proportional to 1 / x dx + J^^ l/xdx + ... = N log [2] , 
where N is the number of decades spanned by the dis- 
tribution (assumed to be a large number). Similarly for 
the digit d, one obtains N\og[{d + l)/d], yielding Ben- 
ford's law. Thus, there is an intimate relation between 
the 1/.T distribution leading to logarithmic, slow relax- 
ations, and Benford's law, stating universal statistics in 
the first digit problem. It remains a challenge to find the 
unifying mechanisms between these various observations. 



THE EXTENT OF APPLICABILITY OF THE 
THEORY 

So far, we have discussed three different physical mech- 
anisms leading to the same aging behavior. There are, of 
course, other forms of slow relaxations in nature, which 
are also commonly observed, and many fascinating phys- 
ical systems which do not fit into the framework outlined 
here. Two important examples are the aging behavior 
of sp in-glas s systems, which appears to be more com- 
plex 5Cll452| and is not described by this model, and the 
behavior of molecular and colloidal glasses, which have 
received much attention in recent years 53|. As indi- 
cated by Refs. [13], [HI], the model proposed here is 
not consistent with the relaxations observed in these sys- 
tems. Another generic form of slow relaxations in na- 
ture is stretched exponential relaxations e^^*/'^-''^, which 
has been experimentally observed in various systems 
[H H, [s^-Ell, and for which several theoretical models 
have been proposed H^-fi^l- For small /3, this form is 
similar to Eq. ([T]), but it is still possible to distinguish 
the two. by analyzing the short-time behavior, where a 
power-law markedly differs from a logarithm. It would 
be interesting to make a better classification of these two 
universality classes (and possibly others), and to deter- 
mine the extent and limitations of the applicability of 
each of them, but this task is beyond the scope of this 
work. 



CONCLUSIONS 

To summarize, in this paper we discussed a generic 
model for slow relaxations and aging, whose signature 
is a distinct crossover from a logarithm to a power-law, 
with no fitting parameters other than the overall scal- 
ing. Theoretically, we have shown how various different 
physical mechanisms give rise to a broad distribution of 
relaxation rates, of a particular form, and analyzed the 
resulting aging behavior. The data from various experi- 
ments was shown to agree well with this prediction, over 
many decades in time. The experiments measure differ- 
ent physical observables, in a variety of systems and tem- 
peratures. This form of aging is fundamentally connected 
to other phenomena which are commonly observed in na- 
ture, such as 1/f noise. Part of the beauty of physics lies 
in the surprising connections it offers, between different 
fields and phenomena. We showed that one could under- 
stand on equal footing the aging of quantum tunneling 
in electron glasses [2^ and the mechanical relaxations of 
plant roots [3], and that both are connected to the 1/f 
noise electrical engineers are well familiar with. As such, 
we believe the model presented can serve as a paradigm 
for slow relaxations and aging for a broad class of sys- 
tems. 
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